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We demonstrate that an appropriate sequence of laser 
pulses allows to condense a gas of trapped bosonic atoms into 
an arbitrary trap level. Such condensation is robust, can be 
achieved in experimentally feasible traps, and may lead to 
multistability and hysteresis phenomena. 



Laser cooling has led to spectacular results in recent 
years So far, however, it has not allowed to reach 
temperatures for which quantum statistics become im- 
portant. In particular, evaporative cooling is used to 
obtain Bose-Einstein condensation of trapped gases . 
Nevertheless, several groups are pursuing the challenging 
goal of condensation via all-optical means . 

In traps of size larger than the inverse wavevector, , 
the temperatures required for condensation are below, or 
of the order of, the photon recoil energy, Er = Tilur = 
Ti 2 k 2 L /2m, where m is the atomic mass. There exist sev- 
eral laser cooling schemes to reach such temperatures 
||[7]]. They exploit single atom "dark states", i.e. states 
which cannot be excited by the laser, but can be pop- 
ulated via spontaneous emission. The main difficulty in 
applying dark state cooling for dense gases is caused by 
light reabsorbtion. Unfortunately, these states are not 
dark with respect to the photons spontaneously emitted 
by other atoms. Thus, at sufficiently high densities, dark 
state cooling may cease to work, since multiple reabsorp- 
tions can increase the system energy by several recoil 
energies per atom p[-ll|. In particular, laser induced 
condensation is feasible only if the reabsorbtion proba- 
bility is smaller than the inverse of the number of energy 
levels accessible via spontaneous emission processes 0j . 

Several remedies to the reabsorption problem have 
been proposed. First, the role of reabsorptions increases 
with the dimensionality. If the reabsorption cross sec- 
tion for trapped atoms is the same as in free space, i.e. 
~ 1/fc 2 , the reabsorptions should not cause any prob- 
lem in ID, have to be carefully considered in 2D, and 
forbid condensation in 3D. Working with quasi-lD or - 
2D systems is thus a possible way to reduce the role of 
reabsorptions The most promising remedy against 
this problem employs the dependence of the reabsorption 
probability for trapped atoms on the fluorescence rate 7, 
which can be easily adjusted in dark state cooling p3| . 
In particular, in the so called Festina Lente limit, when 
7 is much smaller than the trap frequency 10 fl4f , the re- 



absorption processes in which the atoms change energy 
and undergo heating are suppressed. However, neither 
collective cooling schemes in traps of realistic size have 
been investigated in this limit, nor has it been shown that 
laser induced condensation is possible. 

In this Letter we present such investigation. First, 
we formulate the Master Equation (ME) that describes 
Raman cooling in the Festina Lente limit using coarse 
graining in time. Cooling is dynamical, and consists of 
sequences of pairs of laser pulses inducing stimulated and 
spontaneous Raman transitions between the two elec- 
tronic levels of trapped atoms, \g) and |e). The stimu- 
lated absorption pulses induce the energy selective tran- 
sition \g) — > |e) that depopulates all motional states ex- 
cept the dark ones; the spontaneous emission pulses are 
non-selective, and repump the atoms from |e) to \g) pop- 
ulating all accessible motional states. We simulate the 
dynamics generated by the ME in 3D, and show that: i) 
laser induced condensation into an arbitrary trap level is 
possible; an arbitrary trap level may be made dark; ii) 
condensation is robust with respect to changes of physi- 
cal parameters; dark states do not have to be completely 
dark; iii) multistability and hysteresis occur when the 
parameters undergo large changes. In the limit of large 
number of atoms analytic solutions of the ME are found. 

We consider N bosonic atoms with two levels \g) 
and |e) in a non- isotropic trap with the frequencies 
^xyzi^xyz different for the ground and the excited 
states, and non-commensurable one with another. This 
assumption simplifies enormously the dynamics of the 
spontaneous emission processes in the Festina Lente 
limit. We use the coarse graining in time and describe 
variations of the atomic state after one absorption and 
one spontaneous emission pulse. After such cooling cycle 
all atoms are in the ground internal state described by 
the density matrix p(t). This matrix is diagonal in the 
Fock representation corresponding to the bare trap lev- 
els. In order to derive the ME, we separate the effects 
of laser cooling from the ones due to atom-atom colli- 
sions. The latter can be described by a quantum kinetic 
ME, which has been studied in Ref. In this paper 

we concentrate on the laser cooling process. Thus, our 
results are valid in the case when collision processes are 
slow compared to laser cooling. 

For the stimulated Raman transitions, we assume weak 
pulses of duration r a fc s . Their effects can thus be de- 
scribed by second order perturbation theory (formally 
that corresponds to one atom excited at most), 
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where gm-,g\ n ( e ii e j) are bosonic annihilation and cre- 
ation operators of atoms in the ground (excited) in- 
ternal state and in the trap level m = (m x ,m y ,m z ) 
[I = (l x ,l y ,l z )]. The absorption probabilities de- 
scribe transitions from the ground state level m to an 
excited state level /, and depend on Raman laser pulses. 
For instance, if the Raman transition is characterized by 
the maximal (effective) Rabi frequency Oo, temporal en- 
velope /(i), wavevector k^, frequency lul, and detuning 
S = lul — uj a , where uj a is the internal levels energy dif- 
ference, the probabilities are given by 



IY^f = ^|(Z| e ik - x |m}| 2 |/(<5 - wf + OP, 
where the first term is Franck-Condon factor, and 
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f(t)e- ut dt, 



is a function peaked at S = with width of the order 
of 1 /Tats- In the following we use Gaussian pulses with 
f(t) = exp(— t 2 /r 2 bs ) an d UTabs > 1, so that only reso- 
nant Tf£f are relevant (\S - uf + wf„ < t~^ s ). On the 
other hand, fioTahs < 1, so that only a small fraction 
of atoms is excited by the absorption pulse, which en- 
sures the validity of the perturbative approach. In Eq. 
(|l|) two kinds of terms are omitted, since their contri- 
butions vanish after the repumping process: i) the first 
order terms, which correspond to coherences with respect 
to internal levels (off-diagonal matrix elements between 
the states with one excited, and no excited atoms); these 
terms are destroyed after spontaneous emission; ii) coher- 
ences between different excited trap levels; those terms 
are neglected since spontaneous emission in the Festina 
Lente limit is purely diagonal in the I index; thus I ^ V 
coherences are destroyed in the quantum jump down to 
the ground state, and do not influence the dynamics of 
diagonal matrix elements |]l6f . 

The repumping pulses have constant amplitude and 
duration r sp long enough to depopulate totally the ex- 
cited states. In accordance with the Festina Lente limit, 
we will assume that (collective) spontaneous emission 
rates are small in comparison to the trap frequency p| . 
This allows to perform a secular approximation in the 
master equation describing this process, which leads to 
the following equation after spontaneous emission has 
taken place (r ch = T abs + r sp ): 

P(t + T ch ) = p(t) - E 17^ [glg m p(t) + p{t)glg m ] 

I m 
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The last term in Eq. @ describes the integral over all 
possible times r in which the quantum jump from the ex- 
cited state I to the ground state n occurs. The amplitude 
of the excited state is damped during the time r with the 
collective rate ^ fc r^g/effjj, which evidently contain Bose 
enhancement factors, i.e. is proportional to occupation 
numbers plus one of the corresponding ground trap lev- 
els. The integral over r can be extended to 00 since r sp 
is long enough. The spontaneous emission rates are 



r 



sp 



7 / dnW(Q)\(n\e 



ik a n(n)x 



Ol 2 



(4) 



(3) 



where 27 is the single atom (effective) spontaneous emis- 
sion rate (i.e. spontaneous Raman transition rate, con- 
trollable in experiments), integration extends over the 
solid angle f2, W(il) describes dipole radiation pattern, 
n(r2) is a unit vector in the £1 direction, and k a — uj a / c. 

Eq. (||) describes an elementary cooling step, i.e. maps 
a diagonal density operator p(t) onto a diagonal p(t+T c h), 
with all atoms in the ground state. It can thus be simu- 
lated using standard Monte Carlo procedures. We have 
performed such simulations for N = 1 up to N = 500 
atoms, in various dimensions, and for various cooling 
strategies. Franck-Condon factors and trap frequencies 
can be efficiently approximated using the states of an 
isotropic trap of frequency w. We concentrate on 3D 
cooling, beyond the Lamb-Dicke limit, i.e. for the traps 
for which the Lamb-Dicke parameter 77 = ^jEn/huj is 
larger than one. Due to memory storage and calculation 
times, our numerical simulation has to be restricted to 
values of 77 < 2. In that case the atom (having initially 
an energy of the order of few En) may increase its trap 
energy level in the spontaneous emission process by en- 
ergies ~ En, and non-standard cooling schemes have to 
be used to avoid such heating effects O] ■ Generalization 
of the approach of Ref. |l7| allows to cool dynamically 
(i.e. by changing absorption laser pulses appropriately) 
individual atoms to arbitrary trap levels Jl8| . 

The full dynamical cooling cycle must contain se- 
quences of absorption pulses of appropriately chosen fre- 
quencies. The Fourier bandwidth of the pulses can be 
smaller than u>, so that various resonance conditions may 
be employed. We use the following types of pulses: i) 
confinement pulses: spontaneous emission may increase 
each of the quantum numbers m XtVtZ by 0(n 2 ). In D- 
dimensions pulses with detuning 5 = — Df] 2 u), where fj 2 
is the closest integer to n 2 , have thus an overall cool- 
ing effect, and confine the atoms in the energy band 
of D recoils. The use of two sligthly detuned confine- 
ment pulses is recommended; ii) dark-state cooling pulses: 
these pulses should fulfill dark state condition for a se- 
lected state to which the cooling should occur; iii) side- 
band and auxiliary cooling pulses: in general, dark state 
cooling pulses might lead to unexpected trapping in other 
levels. In order to avoid it, auxiliary pulses that empty 
undesired dark states and do not empty the desired dark 
state are needed. For cooling into the ground state, for 
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instance, the sideband cooling pulse with 5 — —uj is used; 
iv) pseudo-confining pulses: with the use of pulses i)-iii) 
cooling is typically very slow. In order to shorten cooling 
time we use pulses with 6 = —3t] 2 uj/2 and S = —t] 2 uj, 
which pseudo-confine the atoms below n — irj 2 /2 and 
n = rj 2 . 
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Population of (0, 0, 0) as a function of number of 
cycles, each consisting of a sequence of absorption pulses with 
s=-12,-6,-4, 0,-13,-7,-5,-1, and repumping pulses. For all the 



FIG. 1. 



absorption pulses A x 



1, except for s — for which 



A z — —2. The initial distribution is thermal with mean 
(n) = 6, and r) = 2.0. Dotted (solid) lines represent the 
case of N = 1 (500) atoms. 



Let us now consider possible dark state conditions, de- 
termined by the total probabilities of emptying a given 
level m, T m — J^i ^tm ■ We cons ider three Raman tran- 
sitions induced by laser pairs propagating in directions 
x, y, and z characterised by three different effective Rabi 
frequencies Qof{t)Aj, where Aj =XViZ account for differ- 
ence of intensities or dephasing between the lasers. As 
we discussed in detail in Ref. Q the dark states may ap- 
pear due to two reasons: vanishing of all Franck-Condon 
coefficients for three lasers, or destructive interference 
effects between the three lasers. The first type of con- 
ditions can be achieved by choosing appropriate detun- 
ing S — slu with integer s: for instance, ground state 
to = (0,0,0) is dark with respect to side-band cooling 
pulses with s = — 1. Since the Franck-Condon coeffi- 
cient (m x + s\e LX \m x ) vanishes for m x — 1 provided 
i] 2 = s + 1, and analogous property holds for y and z, 
the state (1,1,1) can be made dark provided r] 2 is inte- 
ger, and s = rj 2 — 1. Similar, conditions can be found 
for to = (2,2,2). The second type of dark state con- 
ditions correspond to use of resonant absorption with 
s = 0. Chosing for instance A x — 1 and A V:Z such that 
1 + A y + A z = makes all states (to, to, to) dark. 

In Fig. 1 we present our result for ground state cool- 
ing of 1 and 500 sodium atoms in 3D trap with rj = 2, 
using 20 3D-energy shells (i.e. 1771 trap levels). The 
initial state of the system corresponds to mean energy 
QTiu, and is the same for all figures. The pulse sequence 
is s = -12,-6,-4, 0, -13, -7 
s = 0, for which A xv = 1 



A x , 
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1, except for 



-5,-1. 

A z — —2. Pulses 1 and 5 are 
confining, 2,3,6 and 7 pseudo-confining, and 4 and 8 are 
dark state cooling pulses. The many body effects intro- 
duce one very important element to the dynamics: Bose 



enhancement factors, that speed up the dynamics enor- 
mously. The time scale is such that each cooling cycle 
must be, say few times longer than 2-k/uj ~ 10~ 4 s flT9| . 
Then the function / is sufficiently narrow to neglect non- 
resonant transitions. Therefore 1000 cycles correspond to 
about Is. Keeping T fixed, cooling of one atom requires 
here few seconds, whereas collective cooling takes about 
0.1s. After achieving condensation with 500 atoms, con- 
fining pulses can be avoided, a single dark state pulse can 
keep the atoms in the condensed state. 
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FIG. 2. Population of (1, 1, 1) as a function of number of 
cycles, each consisting of a sequence of absorbtion pulses with 
s =-12,-6,-3,3,-13,-7,-4,-2, and A x>VtZ = 1. The initial distri- 
bution is thermal with mean (n) — 6. The cases of r\ = 2.0 
and 2.05, N = 1 and 500 are depicted. 



In Fig. 2 we show results for cooling into the state 
(1, 1, 1) using the sequence s =-12,-6,-3,3,-13,7,-4,-2, with 
A x , v ,z — 1- Here the pulse 4 is the dark state pulse 
(r/ 2 = s + 1), the other are either confining or auxil- 
iary. First, note that when dark state condition is ful- 
filled exactly (77 = 2), cooling of a single atom to (1, 1, 1), 
although slow, is possible. This cooling mechanism is, 
however, very fragile. A tiny perturbation of the dark 
state (77 = 2.05) makes efficient cooling impossible. This 
conclusion does not hold for many atoms though. Quan- 
tum statistics helps to achieve 100% condensation that 
is robust with respect to parameter changes; the results 
for r\ = 2 or 2.05 are undistinguishable, cooling is much 
shorter than in the 1 atom case, and takes about Is pfj[ . 

The existence of various stationary states in our sys- 
tem suggest the possibility of multistability and hystere- 
sis effects [|lj. Indeed, in Fig. 3 we investigate the same 
cooling sequence as in Fig. 1, except that sg=-2, and for 
pulse 4 A z varies adiabatically from —1.94 to —0.08 and 
back during 37200 cycles. For A z close to -2, the ground 
state is in this case dark for s = 0. As A z grows, at some 
point (A z = —2/3) the states (1, 0, 1) and (0, 1, 1) become 
dark. This occurs when the destructive interference 

(OxIe'^'lO.) + (l y \e^y\l y ) = -A z (l z \e lk ^ z \l z ) (5) 

takes place. The system shows multistability and hys- 
teresis: the transfer from (0,0,0) to (1,0, 1) and (0, 1, 1) 
occurs for higher values of A z than vice versa. 
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FIG. 3. Hysteresis cycle in the condensation dynamics into 
(0,0,0), for a sequence of pulses s =-12,-6,-4,0,-13,-7,-5,-2 
with A XlVl z = 1 for all the pulses except for the pulse s = 0. 
For s — A z = —1.94 for the first 1200 pulses, A z increases 
adiabatically during 18600 cycles up to —0.08 (solid line). In 
the next 18600 cycles A z goes back to its original value —1.94 
(dotted line). 

If most atoms are condensed, confining pulses are no 
more needed. One can study then the stationary limit 
of the ME with cooling pulses with fixed s. Amaz- 
ingly, using the Glauber's P representation and sys- 
tematic 1/N expansion the ME can be solved analyti- 
cally. The results can be summarized as follows: i) con- 
ditions for absorption and spontaneous emission rates 
that lead to cooling into an arbitrary state are analyt- 
ically obtained; the system condenses into no-state iff 
f = V (r abs — r abs V sp /V sp ) > O for all n ^ ?? n 

L n — Z^m\ L ran 1 mn rim 1 1 n ra) u iul dn 11 r ll 0- 

ii) the dynamics exhibits two time scales: on a fast scale 
(several cooling cycles), noncondensed modes behave as 
independent quantum harmonic oscillators that approach 
thermal equilibrium (quantum Ornstein-Uhlenbeck pro- 
cesses). The condensate mode is correlated to that dy- 
namics through the atom number conservation. The 
cooling time is of the order of max n (r c ^/r„). On a slow 
time scale (that is N times slower than the fast one) the 
dynamics is dominated by self-transitions from the con- 
densate and back. That produces slow phase diffusion 
and slow decay of the two-time correlation function of 
the condensed mode with the rate ~ J2 m ^mn /^ T ch- 
The scattered photon statistics is Poissonian. 

Sumarizing, using the quantum ME in the Festina 
Lente limit, we have demonstrated that properly de- 
signed sequences of laser pulses allow to condense a gas 
of trapped bosonic atoms into an arbitrary state. Such 
condensation is robust, can be achieved in experimentally 
feasible traps, and leads to multistability and hysteresis 
phenomena. We have neclectcd in our approach atom- 
atom collisions. Collisionally induced population redis- 
tribution should not affect condensate in the (collision- 
ally modified) ground state (0, 0, 0), since this state corre- 
sponds to thermal equilibrium at very low temperatures. 
Thermalisation mechanism might destroy condensates in 
excited states, but that depends on time scales. Conden- 
sation requires seconds, but we have not attempted to 
optimize this time. Optimal times become shorter than 
collisional thermalisation time if N is not too large, and 
T] not too small. 
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